Phase transition in linear sigma model and disoriented chiral condensate 
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We have investigated the phase transition and disoriented chiral condensate domain formation 
in linear sigma model. Solving the equation of motion for the sigma model fields in contact with 
a heat bath, we have shown that the fields undergo phase transition above a certain temperature 
(Tc). It was also shown that when the fields thermalised at temperature above Tc are cooled down 
sufficiently rapidly, disoriented chiral condensate domains are formed quite late in the evolution. 
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- - - The possibility of forming disoriented chiral condensate (DCC) in relativistic heavy ion collisions has generated 

, considerable research activities in recent years. The idea was first proposed by Rajagopal and Wilczek [1~4]. They 
2^ ' argued that for a second order chiral phase transition, the chiral condensate can become temporarily disoriented in 
', the nonequilibrium conditions encountered in heavy ion collisions. As the temperature drops below Tc, the chiral 
, ■ symmetry begins to break by developing domains in which the chiral field is misaligned from its true vacuum value. 
^ ' The misaligned condensate has the same quark content and quantum numbers as do pions and essentially constitute 
; a classical pion field. The system will finally relaxes to the true vacuum and in the process can emit coherent pions. 
QQ ' Since the disoriented domains have well defined isospin orientation, the associated pions can exhibit novel centauro-like 
04 . [5^8] fluctuations of neutral and charged pions [9-12]. 

Most dynamical studies of DCC have been based on the linear sigma model, in which the chiral degrees of freedom 
[ are described by the real 0(4) field <i> — (ct, tt), with the Lagrangian, 
^ . 

c = Ud.'^r - jc^' - f.)\ (1) 

cn : ^4 

, where A is a positive coupling constant and is the pion decay constant. At finite temperature, to leading order in 
A, the thermal fiuctuations < ScjP > of the pions and cr-mesons do generate an effective Hartree type dynamical mass 
\ giving rise to an effective temperature dependent potential AT^/2 [13]. Resulting chiral phase transition is compatible 
■ with the expectations of lattice gauge QCD calculations [14]. One generally introduces a finite symmetry breaking 
' term /lo- to take into account the finite pion mass. However, at present, we are ignoring such terms, as our interest is 
to investigate chiral phase transition and subsequent disoriented chiral condensate domain formation. With symmetry 
O ^' breaking term, there is no exact phase transition. 

1^ ,, Aim of the present letter is to study the influence of external source on the phase transition and subsequent DCC 
i-S^ domain formation when the chiral symmetric phase relaxes back to symmetry broken phase. Indeed, in a heavy ion 
^ collision, while it is possible that in a certain region chiral symmetry is restored, that region must be in contact 
with some environment i.e. background. Exact nature of the environment is difficult to determine but presumably it 
will be consists of mesons and hadron (pions, nucleon etc.). Recognising the uncertainty in the exact nature of the 
environment, we choose to represent it by a white noise source, i.e. a heat bath. To be consistent with fluctuation- 
dissipation theorem, we also include a dissipative term in the equation of motion for the sigma model fields. Thus we 
are studying essentially Langevin equation for the 0(4) fields. Recently it has been shown that in the </>'' theory, hard 
modes can be integrated out on a two loop basis resulting in a Langevin type equations for the soft modes [15,16], 
there by justifying our approach. Langevin equation for 0(4) fields have been used by several authors [17,18] to 
study the interplay of friction and white noise in disoriented chiral condensate formation. We thus propose to study 
following Langevin equation, 

+ (; + -ri^ - £ - ^ - ^ir^I + A(*= - /; - T72)l* ^ C(r,., „, Y) (2) 

where we have used proper time (r) and rapidity (Y), which are the appropriate coordinates for heavy ion scattering. 
In eq. 2 77 is the friction. As told earlier, the environment or the heat bath {(,) was represented as a white noise source 
with zero average and correlation as demanded by fiuctuation-dissipation theorem, 

<ar,x,y,Y)>=0 (3a) 
J < Cairi,xi,yi,Yi)Cb{T2,X2,y2,Y2) > d'^x ^ 2TriSab (3b) 
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Set of partial differential cqs. 2 were solved on a 32'' lattice using a lattice spacing of 1 fm, thus cutting off 
modes with momenta >200 MeV. To show that the model undergoes phase transition we keep the heat bath at 
fixed temperature T and evolve the fields for sufficiently long time (30 fm) such that equilibrium is reached. Fields 
thermalised at higher temperature will be more randomized than fields thermalised at lower temperature. At low 
temperature, the randomisation will not be complete and < a > will have definite non-zero value. Above the critical 
temperature {Tc), randomisation will be complete and < a > will be zero, indicating restoration of symmetry. 

Solving eqn.2 require initial conditions {(j) and (j)). We distribute the initial fields according to a random Gaussian 
with, 

<a>= f{r)U (4a) 

< TTj >= (4b) 

<C7^>-<C7 >2=< nf>-< TT, >2= /2/4/(r) (4c) 

< a >= < TTj >= (4d) 
< &^ >=< n >2= /2/(r) (4e) 



The interpolation function 

/(r) = [f + exp{r - ro)/r)]-i (5) 

separates the central region from the rest of the system. We have taken ro=ll fm and r=0.5 fm. The initial field 
configurations corresponds to quench like condition [1] but it is important to note that the field configuration at 
equilibrium will be independent of initial configiiration. Equilibrium value depend on the heat bath only. We have 
verified this in our code. The other parameter of the model is the friction {rj). In the present paper, we use rj = rj^^ + rja- 
and for r]^ and rjcr use values as calculated by Rischke [19] but once again, its precise value is not of importance here, 
as we are looking for fields at equilibrium. Friction determines the rate of approach to equilibrium. This aspect of 
equilibration was also verified. 

The condensate value of the a field can be considered as the order parameter for the phase transition. If there is a 
phase transition at Tc, the order parameter should exactly vanish for temperatures > Tc, while below Tc it will have 
nonzero values [20]. We calculate the order parameter at the end of the evolution as. 

In fig.l, we have shown the variation of the order parameter with temperature. At low temperature, as expected, 
order parameter is around /^=92 MeV. With increasing temperature, < cr > decreases and become zero around 
Tc=120 MeV, indicating restoration of chiral symmetry. It remain zero at higher temperatures also. The behaviour 
of the fT condensate corresponds to a second order phase transition. The critical temperature Tc=120 MeV is also in 
agreement with mean field calculations. 

Chirally symmetric phase at high temperature will roll back to symmetry broken phase as the system cools and 
temperature drops below the transition temperature. As told earlier, it is has been conjectured [1] that domain like 
structure with definite isospin content may emerge during the roll down period. Numerical simulations of linear sigma 
model indicate that with quench like initial condition, domain like configurations do indeed emerge as the chirally 
symmetric phase roll back to symmetry broken state [21]. However, in heavy ion collision, quench is not a natural 
initial condition. < > and < (j) > are in a configuration appropriate for high temperature but that of < c^^ > 
and < (p^ > are characteristic of a lower temperature. On the contrary, thermalised fields are better suited to mimic 
initial conditions that may arise in heavy ion collisions. Here, < >,<</>> as well as < c^^ > and < <j)^ > are in a 
configuration appropriate for high temperature. 

To see whether domain like structure emerges or not with more appropriate initial condition like the thermalised 
fields, we did a demonstrative calculation. The fields were thermalised at T=200 MeV. We then allow the heat bath 
to cool and follow the evolution of the thermalised fields. Evolution of the thermalised fields now will depend on the 
exact nature of the friction, however we choose to use the same friction as before. We assume the following cooling 
law for the heat bath, 

m = Tol (7) 

with n=l (appropriate for 3d scaling expansion). Assuming that the number density is proportional to the square of 
the amplitude, at each lattice point, we calculate the neutral to total pion ratio according to. 
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Very large or small value of the ratio, over an extended spatial zone will be definite indication of disoriented chiral 
condensate domain formation. At this point we would like to note that contour plot of one single component of the 
pion field as has been shown in ref. [21,22] donot necessarily indicate domain formation. 

In fig. 2 we have shown contour plot of the ratio R in xy plane at rapidity Y=0. The initial random distribution 
(panel a) do not show any domain like structure with very high or low value of the ratio. After the thcrmalisation of 
the fields at T=200 MeV, no domain like structure is evident (panel b). Some small zones with large value of R are 
formed, but they are just thermal fluctuations. As will be shown below, thermal flelds have larger correlation than 
random fields in quenched initial condition. Small zones with large r are manifestation of that correlation. Large 
domain like structure with high/low value of the ratio R starts to emerge after 10-15 fm of evolution and cooling. 
With time domain like structure grow. At other rapidities also, similar behaviour is seen. Fig. 2, clearly demonstrate 
that with appropriate cooling law, multiple disoriented chiral condensate domains with different isospin orientations 
can form as the chirally symmetric phase roll back to broken phase. However domain formation occur quite late in 
the evolution. With the scaling cooling law, by the time domain like structures emerge, the system is cooled to ~20 
MeV. It is doubtful whether in heavy ion collision, system can be allowed to be cooled to this extent. Current wisdom 
is that the hadrons freeze-out around 100-160 MeV. With this reservation in mind, it may be said that even with 
thermal fields domains of disoriented chiral condensate can form. 

Above behaviour is confirmed from the correlation study also. We define a correlation function at rapidity Y as 
[21], 

where the sum is taken over those grid points i and j such that the distance between i and j is r. In fig. 3, we have 
shown the evolution of the correlation function at rapidity Y=0. Initially there is no correlation length beyond the 
lattice spacing of 1 fm. After thcrmalisation, the correlation is increased marginally. At the thcrmalisation, the fields 
arc at the minimum of the potential, consequently develops some correlation. Correlation donot increase till 10 fm of 
evolution. At later time, long range correlation develops, increasing with time. At 20 and 30 fm, pions separated by 
large distance are correlated. 

If a single DCC domain is formed in heavy ion collision, it can be easily detected. Probability to obtain a particular 
fraction f of (say) neutral pion from a single domain is P{f) = l/2y/J [1]. However, our simulation indicate that a 
few numbers of domains with different isospin orientations are formed. Naturally the resultant distribution will not 
be l/^/f type. In fig. 4, / distribution for 100 events, at zero rapidity are shown. Density of pions at rapidity Y was 
calculated by integrating over the space-time as, 

N^{Y) = j TT^TdTdxdy (10) 

If the thermalised fields are evolved upto r=10 fm, then the / distribution is sharply peaked around the isospin 
symmetric average value of 1/3. At early stage of evolution, as indicated in fig. 2 and 3, long range correlations or 
domain like structure are not developed in individual events. If the fields are evolved upto 20 or 30 fm, then there is 
definite domain like formation and also long range correlation is developed, in individual events. The / distribution 
still has average /13, but the width is increased. This indicate that large fiuctuations in the ratio r from event to 
event is expected when there is DCC like phenomena. Experimentally, those events for which neutral to total pion 
ratio shows considerable fluctuations will be promising candidate for DCC events. 

In summary, we have studied the phase transition and disoriented chiral condensate formation in linear sigma model 
in presence of background. The background was represented by a white noise source i.e. a heatbath. Evolving the 
fields in contact with the heat bath for sufficiently long time such that fields are thermalised, it was shown the the 
model undergoes 2nd order phase transition around Tc=120 MeV. It was also shown that if the thermalised fields are 
cooled down sufficiently rapidly, multiple disoriented chiral condensate domains are formed. Long range correlation 
also develops. However, even with sufficient fast cooling, with thermalised fields, domains or long range correlations 
develop quite late in the evolution (after 10-15 fm of cooling, when the fields are cooled to 20 MeV or so). It is 
doubtful whether in heavy ion collision, the hadronic fluid can be allowed to evolve and cool for such long time. It is 
possible that they undergoes freeze-out much earlier. Then, disoriented chiral condensate domains will not be formed. 
We have also shown that with multiple domain formation, probability distribution donot follow the l/\/J type of 
law, rather, it is a Gaussian with average at the isospin symmetric value of 1/3, but with increased fluctuations. It 
is suggested that one should look into events with large fluctuations of the neutral to pion ratio around the isospin 
symmetric value of 1/3. 
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FIG. 1. Equilibrium value of the order parameter as a function of the temperature. 




FIG. 2. contour plot of the neutral to total pion ratio at rapidity Y=0 at different stages of evolution. 
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FIG. 4. Probabilty distribution for the neutral to total pion ratio at rapidity Y=0, for 100 events at different stages of 
evolution 
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